RENORMGROUP INVARIANTS AND APPROXIMATIONS OF MAPPINGS 



Gennady N. Nikolaev 
Institute of Automation and Electrometry of SB RAS, 
Pr. Koptyuga 1, Novosibirsk, 630090, Russia 
(Dated: November 18, 2004) 

The relationship between mappings ol sets and renormahzation group transformations is estab- 
lished, and renormalization group invariants of such mappings are found. These results are valid 
both for continuous and discrete mappings and for various dimensionality of image and preimage 
of the mappings too. It is suggested a number of mapping approximations improved in comparison 
with an ordinary power expansion. The approximations take into account the global one-to-one 
character of the mappings. The method is illustrated by a number of examples: by reconstructing 
' of some analytical functions, calculating the integral of the typical partition function of statistical 

, mechanics and the ground state energy for the quartic anharmonic oscillator. In the whole range 

• of nonlinearity parameter varying from zero up to infinity the accuracy of the RG approximation 

04 ' based on a few terms of divergent series is about 0.06% in the next to last case and 0.004% in the 

, last one. 

Q . PACS numbers: 02.70.-c; 02.30.-f; 02.30. Mv; 05.20.-y; 01.90.+g 

Q> . I. INTRODUCTION 

Renormalization group (RG) conception originally arose in the beginning of the fifties I*! as a result of the discovery 
^ , of special group of continuous transformations in quantum electrodynamics. This transformations are connected with 
ly-j ■ a complicated procedure of renormalization, that is 'removing of ultra-violet infinities'. It turned out, that physical 
\^ [ quantities do not vary at simultaneous rescaling of the 4- momentum transfer squared z, mass squared y, and also 
> special transformation of the dimensionless charge squared g {g ^ 'g{z, y, g)). The invariant charge 'g{z, y, g) obeys to 
the equation 0,0 

9{z,y,g) = g{z/t,y/t,g{t,y,g)). (1) 
For the massless quantum-field model instead of ^ the following relation is valid 

g{z,g)^g{z/t,g{t,g)). (2) 
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The generalization of equation Q for the two-charge quantum-field model was proposed also In this case, instead 
of expression Q the following two relations are fulfilled 

9i(.z,gi.,g2) ^gi(zlt,g^{t,gi,g2),g2{t, 91,92)) i ^ 1,2 (3) 

The general solution of such functional equations was obtained 0] . 

The rapid expansion of the RG method far beyond the quantum-field theories began since the seventies. It was 
strongly promoted by the article Q . In this article the conception of an approximate RG was introduced for analysis 
of thermodynamic systems behaviour near the phase transition points. Afterwards the RG method was applied for 
description of turbulence |^, , polymeric compounds , radiation transfer in opaque mediums with the strong 
frequency dependence of the quantum path length Q , fractals [1(tI | , scripts of the dynamic systems transitions to the 
determined chaos [TH . and also was applied to other problem. The relation was found out between nonlinear problems 
of radiation transport and additive variant of the RG |0, . This result was extended on the wide range of the 
physical problems; and the functional self-similarity concept, that generalizes the usual self-similarity concept, was 
introduced 0|. The modern state of the researches concerned with use of this concept for boundary- value problems 
of the mathematical physics is represented in the recent review 

In the present paper a relationship between mappings from one set into another one and RG transformations is 
established, and RG invariants of such mappings are found. On this basis a number of RG approximations that are 
founded on a few terms of series of one-to-one functions is constructed. They represent better approximations of the 
sought function in comparison with its Taylor series expansion because the RG-approximations inherits its one-to-one 
nature. 

The necessity for reasonable approximations of physical quantities is connected to the fact that the majority of 
realistic problems cannot be solved exactly. To solve them, one uses various approximating methods that are in general 
different kinds of perturbation theories near a chosen zero approximation. These methods are evidently convenient if 
the perturbation (interaction or coupling) parameter g is small, (7 <C 1. However, if 5 ^ 1 the expansion in powers 



2 



of g have no sense. But for many realistic systems, because of their complexity, it is technically impossible to use 
other calculation techniques except for the perturbation one. In such a case there are a number of resummation 
(or reconstruction) methods that allow to find the value of some function f{g) in the case oi g ^ 1 using only the 
results of its calculation in the limit g <^ 1. These methods include the improved perturbation theory, the Pade 
approximation, the Borel summation, the conformal mapping and their combinations |l6lll7| . the continued fraction 
approximations and the Romberg algorithm 19|. These methods, to be accurate enough, need the calculation of 
many subsequent approximations for the sought function. But if system is so complex that we are able to calculate 
solely a couple of approximations, then the majority of these methods lose their applicability. In this situation the 
method of RG approximation seems especially relevant because of its quite good accuracy when using only a few 
terms of the perturbation theory. 



II. RENORMGROUP INVARIANTS OF ONE-TO-ONE MAPPINGS 



At first let us consider a function of a single variable f{x) as a primary mapping. It is assumed hereinafter that 
the function f{x) is one-to-one, i.e. it has the inverse function. Now we consider the value of this function at the 
argument xi + x. Let us express the quantity x via the inverse function: x = f^^{f{x)). This relation allows one to 
write down f{xi + x) as follows: 

fix, +x)^ fix, + f-\fix))) ^ Fix,; fix)). (4) 

Thus, we have presented fix, + x) as a new function F of two variables, value of the initial function / at the original 
argument x and the translation x,. According to definition Q), the value of the function F a,t x, = coincides with 
the value of the initial function / at the argument x: 

FiO-Jix)) = fix) (5) 

Now let us consider the value of the function / at the argument x = x + x, + X2 displaced sequentially on x, and 
X2 with respect to x. The value of the function fix) can be represented via the function F by two equivalent ways. 
According to the one way, it is possible to consider the two sequential displacement x, and X2 as a whole displacement 
Xw = xi + X2 and to take advantage of the definition that leads to the relation 

fix+ixi+X2)) = Fix-Jixi+X2)). (6) 

And then it is possible to take advantage of definition |0J once more to transform the function fix + X2) on the 
right-hand side of expression © in terms of F: 

Fix;fixi+X2))=Fix;Fixi-Jix2))). (7) 

On the other way, at first one can use definition Q with replacement x + xi for x, that results in 

fiix + xi)+X2)^Fix + xi;f{x2)). (8) 

Equating the right parts of expressions Q and ||SJ|, taking into account equality lO, we get the functional equation 

Fix + xi-J))^Fix;Fix,-J)). (9) 

Hereinafter we have omitted for brevity the argument x of the function /. Equation (^J coincides in its structure with 
the equation ^ for the radiation transport in [T^ . 

Note that relation @ turns into Q by means of the replacements x = ln(z/i), xi = ln(i), and the obvious change 
of the notations: f ^ g, F ~^'g. 

From equation the partial differential equation for the function F is easily obtained For this purpose 

it is enough to differentiate both parts of equation Q on xi and to take advantage of property ||SJ): 

[-^-Pif)-^]Fix;f)^RFix-J)^0; (10) 
dFixiJ) 

Pif) = Ui=o- (11) 

Functional equation (|5J) expresses the invariance of F under the RG transformation |l4l Il5j representing a simulta- 
neous one-parametric point transformation of both arguments of the function F: 

Tixi): x^S: = ~xi+x, f ^ f ^ Fix^J); F(0; /) = /. (12) 

Equation @ guarantees fulfillment of the group property T(2;i)r(x2) = Tixi + X2)- Note that differential equation 
(|10|l represents the infinitesimal form of transformation \V2.\ . and R is the infinitesimal operator of RG-symmetry. 
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III. RENORMGROUP APPROXIMATIONS OF ONE-TO-ONE MAPPINGS 



Frequently it is possible to define the behavior of some sought function only in a vicinity of some point. Such 
situation is typical for iterative procedures or for one of the basic methods of theoretical and mathematical physics, 
perturbation theory. As a rule, the calculation of the higher-order corrections is connected with significant difhculties 
of the technical or basic character (for example, because of the divergence of the higher-order corrections). It is 
known also that a finite series function approximation often describes wrongly the long range function behaviour. 
Relations ©, l|lt)|l allow to receive improved approximation of the sought function (so-called RG method This 
approximation more correctly describes the long range behavior of the sought functions. The solution of equation 
(|10|l with boundary condition ^ is possible to represent as 

F{x;f) = X-Ux + X(f)), (13) 



W)^/^, (14) 



where is inverse to X function. Ordinarily, to determine the function (|14() one use an approximate value for /? 

calculated in the small parameter limit / = . Let us consider other approach that does not presuppose the 
smallness. According to definitions and ^ we have: 

/3(/)^f ^/(x). (15) 
ax 

Suppose that we know an approximate expression of f{x) in a neighborhood of some point xs, for example, as a 
Taylor series (without loosing of generality one can set xq = 0): 

f{x) = ao + aix + a2X^ + 03^^ + . . . . (16) 

Series can be inverted: 

x^(f> + + h:if + . . . , (17) 

where 4> — {f ~ ao)/ai, the coefficients bi are expressed via the coefficients of initial series (jlSf) . The sequential 
substitution of expansion H16() and (|17|l into the expression (|15(l enables us to get the required dependence of (3{f) as 
an infinite series on </). If we use various finite series approximations of (3{f) , we receive various RG approximations 
of the sought function by the inversion of expression (|14|l . So, the linear approach of (3{f) gives the following RG 
approximation of the sought function f{x): 

2 

fix) ~ fj^\x) = ao -f ^ [exp(2a2a;/ai) - 1] . (18) 
Za2 



From the square-law approach of /?(/) another RG-approximation follows: 

^ \ (19) 



/(x)^/f (a;)=ai + ^ 
' 02 



— KCOth(Kx) 

a2 



where k = \/3[(a2/ai)^ — as/ai]. 

It is interesting to note, that expression H14I) represents exactly the variable a: as a function of / taking into account 
(|15|1 . Thus, expressions (|14|l and IjlTfl are equivalent. Therefore, one can get one more variant of the RG function / 
approximation, using only a few sequential terms of series H17|l and inverting the obtained approximate equation with 
respect to /. So, the square-law approach of series H17|l results in the following RG approximation: 

fix) ^ f^^\x) ^ ao + ^[1 - Vl-4(a2/ai)x] = ao + , ^"^^'^ ^ ^ ■ (20) 

2a2 1 + ^1 - 4(a2/ai)x 

It is possible also to receive other variants of the function RG approximations, using other approaches for fix), 
f ix), or xif) instead of the Taylor series . For example, one can get another RG approximation of the sought function 
using second convergent of the continued fraction corresponding to series 1)17(1 : 

/(.)-/^i/x)=ao + ^3f^ (21) 
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It should be noted that this expression coincides with the approximate one obtained from 1)20(1 as a result of the square 
root expansion. It is interesting to note also the coincidence of this expression with the self-similar approximation 

Let us represent the estimations of the RG approximations accuracy of some functions for x = 1. We use the 
notation Aa — [/a(l) — x 100%, where fa{x) is any of the RG approximations of a function f{x), as a 

(2) (3) (2) (2) 

measure of accuracy. For the function f{x) = In (1 + x) we have « —8.8, Aj^ « 1.6, A)^ fa 5.3, A\-^j « —3.8, 

(2) (3) 

A)p w —27.8, A^^ w 20.2, where the last two quantities refer to the finite Taylor series approximations of the 
square-law type and cubic one respectively. The accuracy of the different RG approximations for other function 
f{x) = x/VTTx is represented as follows: A^^^ « -10.3, Af^ « 3.4, A^^ w 3.5, A^^^ « -5.4, A^^^ « -28.9, 
A^^^ « 24.0. 

It should be noted that the function RG approximations for ai = can be obtained in the same way. In this case 
x{f ) is represented as a power series on 0i = \/ {f — o-o)/a2- 

It is important to note, that the expansion in powers of x of both sought function f{x) and its RG-approximation 

fa"\x) coincide up to the order n inclusive. Therefore a successive RG-approximations f^\x) of the sought function 
f{x) converge to the latter one not worse of the Taylor series while value n sequentially increases. RG-approximations 
of one-to-one function represent in general its better approximations in comparison with its Taylor series expansion 
because RG-approximations inherits its one-to-one nature. 



IV. PARTITION FUNCTIONS OF STATISTICAL MECHANICS 



Let us illustrate accuracy of our method on the two important physical examples. Both of them satisfy the following 
conditions: (a) the standard perturbation theory is not applicable at all leading to divergent series; (b) some analytical 
expressions are known; (c) analytical or numerical calculations are available giving the possibility to compare them 
with obtained results; and (d) there are other approximations that may also be compared with the method considered. 

First example is an integral having the mathematical structure typical for the partition functions for statistical 
problems with effective potential of the form V{4)) = cfP' + gcj)"*", that is, for the problems with the so-called (p^ 
interaction. Consider the integral 



1 r 

l{9) = -r / 



exp(-(/)2 - a<i>'^) d(j) 



with g e (0,oo). The generating functional of the 4)^ quantum field theory has also the structure of H22() . 
Integral H22|l admits an exact quadrature expressed via the modified Bessel function of the second kind: 

exp(l/8ff) 1 
It is easy to find the expansion of in powers of g: 

n») = £tgr(2„4),". 

n=0 V \ / 

The strong coupling limit g 1 of H22|) can be obtain by using saddle-point method 



I{9) 



■2^ 



r(i)..-/^-r(i),-3/4 



> 1 



Substituting here the quantities F (1/4) « 3.624 409, F (3/4) sa 1.225417 for gamma- function F(y), we get 

I{g) ~ 1.022 735 g'^'^ - 0.345 384 g-^'^ ; g > 1 . 
As can see, the nth term in the expression H24I) diverges at any non-zero 17, as n — > 00, since 

F(27i+l/2) /4n\" 



(22) 



n — > c» . 



(23) 



(24) 



(25) 



(26) 



(27) 
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So, expansion in powers of g has no sense aX g ^ 1. Retaining three terms of series (|24() . we get the approximation of 

(|12J) 

J(5)^l-|g + ^5^ (28) 

The accuracy A^-* of the approximation H28|) is about 360 % at 5 = 1. On the other hand, we have A^^-* w —18%, 

A)^ w 7%, and A)^^^ w 9% for RG approximations p8|) . (|20|) . and (|21|l respectively that are based on approximation 
(|28|l . It is seems satisfactory accuracy for approximations based on a few terms of divergent series. 



A. Modified RG approximation. 



One can get more precise approximation by our method combined with the principle of minimal sensitivity (PMS) 
PH even in the whole range g e (0, 00). 

Let us rewrite partition integral H22I) in the equivalent form 



1 r°° 

I{q) = / CXp(-z2^2 ^ ^y) 

Vl" J-oc 



where 

and z is trial parameter. Expanding the integrand of (|29|l in powers of (|30|l . we get 



1 



h{g,z) = /i(g,2) + - 



i5_ 

1 



15.g 

8^ 



10532 
"32^ 



in the first and in the second order respectively. 

Let us make the substitution a; = 1 — l/z^. Then expressions H31II and (|32() transform into 



Ii{g,x) = ^/l 



h{g,x) = Ii{g,x) + \/l-x 



13, ,0 
1+2^-45(1-2;) 



155 



105 
"32" 



Exact function 1(g) does not depend on z (or x) according ly to (123) or (EHIl, therefore ^iM 



(29) 

(30) 

(31) 
(32) 



(33) 
(34) 



g„ = 0. However, any finite 

series approximation In{g,x) actually explicitly depends on x. According to the PMS ,21,], let us require minimal 
sensitivity to x as early as in the first order: 



dli{g,x) 
dx 



that gives the relation 



2 X 



" 5 [l-xf ■ 

We see from here that < x < 1 owing io Q < g < 00. Inverting this equation results in 

1 + 55- VI + IO5 



Substituting ^ into we obtain 



x{9) 



I2 {x) = Vl - a; 



55 



■ 1 3 

^ + r+20" 



(35) 



(36) 



(37) 



(38) 



6 



Finally, after substituting RG approximation (|18l) instead of the square brackets in (|38|l . we get 



lf\g)^./T^W)\l + 



15 



exp ( -^^x (g) 



(39) 



where x{g) is given by formula H37|) . In accordance with the general property of the RG approximation method, the 
weak coupling limit (<? <C 1) of (|39|l up to the second order coincides with expansion (|28ll . The strong coupling limit 
of (123), when 5 > 1, is 



I^ig) ~ 1.164 456 g-i/''- .634 951 g"^/''; g > 1 



(40) 



The accuracy of the approximation H39|) is the monotone increasing function of g that tends to limit « 13.8% when 
00. 



B. Improved modified RG approximation. 



It is important to note, that relation (|36|l was obtained by applying the PMS to the first order approximation l|33() 
but not to the RG approximation. Therefore the relation H36|l is not optimum for the RG approximation to achieve 
the best approximation of the exact function I{g). Let us use more general relation instead of l|36(l 

(41) 



5p (1 — x)^ 

with trial parameter p. Inverting this equation gives 

l + 5pg- VTTTUpg 



5pg 



Using the relation H41I) results in the following RG approximation 



lP{g) = VT^^){l + 



2 [1 + |(P-1)]^ 
15 1 + lip- if 



(42) 



(43) 



where x{g) is given by formula 142|) . 

The weak coupling limit of (|43|l again coincides with the expansion (|28|l up to the second order. As can be shown, 

this remarkable property is the consequence the fact that the modified perturbation potential AV H3U|) approaches 

the initial perturbation potential gcj)'^ when x —> 0. 

On the other hand, the strong coupling limit of H43(l depends on p. The optimal value of the parameter p can be 

(2) 

found now from the coincidence of the strong coupling limit of (g) with this limit of the exact integral I{g) that 
is given by ((JSJ, i.e. 



lim 



I? (9) 



= 1 = 



1/4 



1 



2 [i + i(p-i)] 



2 r 



exp 



3 1 



(P-If 



2p[l + |(p-l 



^(5) ^ \5pJ ' 15 l + |(p-l)2 

Solving the last equality of (|44|l with respect to p we get 

p w 1.779 643. 

Using this value, we get the following strong coupling limit of the RG approximation 143() 



1.022 765. 



r^' (g) ~ 1.022 765 g^^/"* - .343 514 g~ 



3/4 



5> 1- 



(44) 



(45) 



(46) 



By numerical calculations, the maximal error for the RG approximation (|43|l with parameter p given by H45|) is 
found to be « 0.06% at 5 = 1.9. Remind that when we did not optimized p {p = 1), the accuracy of the method was 
13.8%. So, by using asymptotic constraint condition (|44fl we improved the accuracy more than in two hundred times. 
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V. QUARTIC ANHARMONIC OSCILLATOR 



Now let us check the abihty of our method to approximate the ground state energy Eq of the one-dimensional 
anharmonic oscillator with the Hamiltonian (h = I) 



H 



2m dy^ 



(47) 



in which m, lu and A are positive constants, y G (— oo , oo). 

Actually, a lot of various physical models can be reduced to the Hamiltonian (|47|l . The close connection of the 
quartic-oscillator model with the so-called t/)^ model in the quantum field theory is also well acknowledged. 

The Hamiltonian 147(1 , making the change y —f ^/my, can be transformed to the form 



1^ 
2dy^ 



(48) 



which is more convenient for calculations. The perturbation series with respect to the coupling constant A diverges 
for any finite value of this constant [2^ |23| . The divergence of the series for the dimensionless ground state energy 



e(.g) EE Eo/uj {g = X/uj^) 



(49) 



is so strong that the expansion in p owers of g has no sense at g ~ 1. This is due to the increase of coefficients near 
5" according to S^n! as n ^ 00 f23!|. The expansion of (|49|l in the weak-coupling limit USES is 



e(3) 



1 



21 



333 



30885 



79 



-9 



-9 



2 4" 8 16 128 

The asymptotic behavior of (|49|l in the strong-coupling limit is given by HE 1^ 

e(5) fo.667 986 + 0.143 675"^/^ - 0.008 85-^/^ 



(5 « 1) . 



(9 » 1) 



(50) 



(51) 



To estimate our RG approximation, we take an advantage of exact numerical calculations |25[l26il27| for the ground 
state energy (|49|l for a wide range of coupling constants. Retaining three terms of series (|5U|I , we get the approximation 
of (|49|l . the relative accuracy of which is about 270% ai g — I. In turn, RG approximations (|18|) . (|20|l . and l|21|) 



that based on this approximation give us A 



(2) 



25%, A 



reasonable accuracy for approximations based only on two first terms of divergent series. 



—0.5% and Aj^-^y « 17% respectively. This is rather 



A. Modified RG approximation. 

Now we again improve our method to be applicable to arbitrary strong coupling. For doing this we take as a trial 
approximation the Hamiltonian 



1 d2 LJo^ 2 

''' = -2d^''—y ■ 



(52) 



and calculate the ground state energy for 1(491) using the perturbation theory over H — Hq = (o;^ — a;o^)y^/2 + Ay^. 
In the zeroth order approximation we have eo(g,a;o) = ujq/2u}. The first approximation for ((49(1 is 



ei(ff,a;o) = ^ 



Wo/ \UJqJ 



(53) 



and the second one is 



e2{g,uJo) = ei{g) - 



UJo_ 

Wlu 



6.9 



LUo 



(54) 



Let us make the substitution x = 1 — (uj/ujof . Then expressions (|53|l and H54() become 



[l-x) 



-1/2 



(55) 



16 



x-^^il-xf^g +&{l-xfg^\{l~x) 



"1/2 



Applying again the PMS to the first order approximation 

dei{g,x) 



dx 



= 0, 



we get the following relation 



5=7 



6 (l-x)3/2 ■ 

We see from here that < x < 1 because of < g < 00. Substituting H58(l into (|56(l we get 



62(2;) = ^ 



1 



1 



1 — tX — —x'^ 



{l-x) 



-1/2 



4 48 

At last, replacing the square brackets in H59() by the RG approximation (|18|l . we get 



1 r 3 

2 1-2 



exp ( (g) ) - 1 



[l-:^(5)]"'/'- 



Here a;(g) is positive and continuous solution of the equation H58|) with respect to x: 



xig) 



l-|(5c/5)' [2cos(a/3)-l]% 
l-|(5c/<?)' [^+ + ^--1]', 



9 <9c 
9>9c 



where 



A, 



9^/3' 

arccos ^2 [g/gcf ~ l) , 

-l + 2{g/g,f±2{g/g,) \{g/g,f-l 



1/2 



1/3 



(56) 



(57) 



(58) 



(59) 



(60) 



(61) 



(62) 



In accordance with the method of the RG approximation, the week coupling limit of H60() coincides with the 
expansion (|50|l up to the second order . The strong coupling limit of (|60|) . when g 3> 1, is 



e^^^ (5) ~ 5^/3 (.661395 + .148 0353"^/^- .010 255g~'*/^j ; (g > 1). (63) 
The accuracy of the approximation (|60|l monotone decreases with g and approaches a limit ~ 0.99% when g — > 00. 



B. Improved modified RG approximation. 

To improve the accuracy we use again more general relation than (|58|l 
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1 



Qp (1 - a;)3/2 

Using this relation results in the following RG approximation 



2 1 + 6(p-l)2 



ll + 6(p-l)2 
^"P'6 p(2p^l) "(g) 



[1-^(5)] 



'1/2 



(64) 



(65) 
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where x{g) is given by the same formula H61|) as before but now gc = l/(9-\/3p) instead of H62() . The weak coupHng 
hmit of Hfi5(l coincide again with the expansion (|5()|l up to the second order as in the previous section Hvl and on the 
same reasons. On the other hand, the strong couphng hmit of H65(l depends on p. The optimal value of the parameter 

(2) 

p can be found now from the coincidence of the strong coupling limit of ej^ (g) with known exact limit of e{g) which 
is given by (|51|l : 



=,(2) 



lim 



(g) 



= 1 = 



1/3 



1 - 



3 (2p-l) 



e{g) ^ \A J 2 1 + 6(p-l)2 

Solving the last equality of (|66|l with respect to p we get 



exp 



ll + 6(p-l) 
6 p(2p-l) 



- 1 



0.667 986. 



p w 1.472 032. 

With this value we get the following strong coupling limit of the RG approximation 1)65(1 

e{g) ~ (o.667 986 + 0.143 62.g-2/3 _ 0.0O8 65"''/^) ; (ff > 1) 



(66) 



(67) 



(68) 



By numerical calculations, the maximal relative error for the RG approximation (|65() with parameter p given by 
(|H7| is found to be A^,^' « 0.003 6% at g = 0.7. Remind that whe n we did not optimized p (p — 1), the accuracy of 
the method was 0.99%. Thus, we improved the accuracy in about three hundred times using asymptotic constraint 
condition ((^ . 



VI. CONCLUSION 



The method of RG approximation presented in this paper possesses the following peculiarities. It seems to be 
logical and simple. It needs minimal information, e. g., only two first terms of perturbation theory. Even in the 
simplest form the error of its approximation is about 10% at g = 1 even for a sought function that has divergent 
series in powers of g. The accuracy of the modified RG approximation that uses the PMS to make it applicable for 
arbitrary values of g is practically the same. The the accuracy of the approximation can be drastically improved 
in hundreds times while using the PMS if we take advantage of some fitting parameter to be defined by additional 
constraint condition, e.g., by making equal meaning of the sought function and its approximation at g 00. So, the 
accuracy of the RG approximation of partition functions integral H22(l is improved in more than two hundred times 
and is of 0.06%. The accuracy of the RG approximation of the ground state energy of the one-dimensional quartic 
anharmonic oscillator is also improved in the same way in almost three hundred times and is of 0.003 6%. 

The accuracy of our method can be compared with that of other anal ytic al approaches that use an equivalent 
number of approximate terms, that is with the methods |2ft I21L I2M l29l l3(l l3]| . The first two methods are called the 
modified perturbation theory, and the others are called the self-similar app roximation. All of them also use either the 
PMS condition similar to l|35|) , H57I) or principle of minimal difference |2(]| to have an opportunity to approximate an 
interesting physical quantity in the whole range (g G 0,(xi) of parameter g. 

The accuracy of approximation of partition functions integral (|22|l by means of both the self-similar approximation 
[30| and the modified perturbation theory H38|) is about 10%, practically the same as in our method. If in the latter 
one we use instead of (|36ll generalized relation (|41|l with the fitting parameter p being defined by asymptotic constraint 
condition similar to (|44|) then the accuracy of this approximation is improved up to 0.19%. As to the self-similar 
approximation, its accuracy is improved up to 1% by use of the "fixed-point distance" fitting parameter being also 
defined by asymptotic constraint condition. Remember that the method of improved modified RG approximation 
presented in this paper has the accuracy 0.06%. 

In the case of calculation of the ground state energy of the one-dimensional anharmonic quartic oscillator the 
accuracy of the self-similar approximation is about 10% while using the principle of minimal difference and is 0.3% 
if using the PMS '29','3T| (the accuracy cited in '2^ corresponds to only strong-coupling limit). Different variants |^ 
I2S.J -32,1 ii2o (see also formula (|59|l 'l of the second-order modified perturbation theory give the ground-state energy 

with the accuracy not higher than 1%. Notice that if we substitute in the second-order approximation H56() ') instead 
of (|58|l the generalized relation H64|l with the fitting parameter p being defined by asymptotic constraint condition 
similar to H66I) then the accuracy of this approximation is improved up to 0.008%. Remember that the accuracy of 
our modified RG approximation (|60|l is about 1% whereas the accuracy of improved modified RG approximation l|65|) 
is 0.003 6%. 

Our approach allows a number of generalizations. 



10 



In general, the invariance relation (O with respect to the RG transformation H12(l is valid for any non-degenerate 
mapping of an associative set of elements x in other set. Indeed, relation @ was obtained with the assumption that 
the variable value set a; is a semigroup with respect to '+' operation. Evidently this set can be either continuous or 
discrete. In the latter case it is possible to represent the successive approximations of any function / as some discrete 
mapping. Let us use the number of the approach (0,1,2,...) as the argument of the discrete mapping and the 
successive approximation of the fmiction (/o(2:), /i (a;), f2{x), . . .) as the value of the mapping. In this case invariance 
relation Q reads: 

F{n + p- fo{x)) = F{n; F{p; fo{x))), (69) 

where n and p are non-negative integers. This relation 1()9|) can be rewritten in the other form, if one use these discrete 
variables as the subscripts of the mapping F: 

F,,+piMx)) ^ F^iFpiMx))). (70) 

Such functional self-similar transformation was obtained pol | from other reasons and also was used for the constructing 
of the so called self-similar approximation. 

Generally a mapping F connects elements x of one set M of the dimension m and elements / of another set N of the 
dimension n. In addition the RG invariance relation @ is valid in the case of one-to-one correspondence between the 
image and preimage of the mapping. For the continuous sets and differentiable mappings this requirement is fulfilled, 
if the rank of the functional determinant 9(/i, /2, . . . , fn)/d{xi, X2, ■ ■ ■ , Xm) is equal to the least of integers m and n 
(so called nondegenerate mapping). In the case of m < n the nondegenerate mapping F defines some m— parametric 
nonintersecting subset in the set N (some m— dimensional nonintersecting hypersurface) . The particular case m = 1 
and n = 2 corresponds to the two-charge quantum- field model . In this model the invariant charges obeys functional 
equations In case of m > n the mapping definition domain that is invariant with respect to RG transformation 
can be anyone n— dimensional nonintersecting hypersurface in M. 

One should notice that the well known monotonous dependence of the invariant charge on the 4-momentum transfer 
squared in the gauge theories of the electro- weak and strong interactions that are renormalizable seems natural in the 
light stated above. 

In this paper the RG approximation presented is based on the expansion of some function f{x) in natural number 
powers of (x — xq) around some point xq. However, our approach can be generalized to the case of both meromorphic 
functions that allow their expansion in positive and negative integer powers and arbitrary (nonmonotonic) functions. 
The latter generalization consists in representing of a real function as the real or imaginary part of some complex 
one-to-one function and RG approximating of the latter one. 

These and others generalizations of our approach will be given in further publications. 

The author would like to gratefully acknowledge the useful discussion with S.G. Rautian. 
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